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ON THE ACTUAL AND PROBABLE ERRORS OF INTERPO- 
LATED VALUES DERIVED FROM NUMERICAL TABLES BY 
MEANS OF FIRST DIFFERENCES. 



BY E. S. WOODWAED, C. E. 
[Continued from page 149.] 
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(a — b) will be the probable error or r when the relation between a and b is 
such that 

jr- = J, or when b = %a. 

In this case r = a — b = a — J<x = \a. When b > \a and < a, r will 
be given by (a) ; and when b < Ja, r will be given by (/?). 

Evaluating (a) there results r = ± [a -f- 6 — j/(2a&)], while (/?) gives 
r = ± Jo. Therefore, when 6 lies between and |a 

r = ± Ja; (11) 

and when 6 lies between \a and a 

r = ± [a + b — i/(2a5)]. (12) 

In the problem of interpolation, a = (1 — 1)\ and b = t.\. Supposing 
l-t>t, 

Therefore, according to equations (11) and (12), 

when t lies between and J, r = ± £(1 — t), ~) 

" t" " $ « |, r = ± l[l-i/2p=T)], V (13) 

" < " " f " 1, r = ± Jfc J 

It appears from these formulas that the probable error of an interpolated 
value of the first species is always less than the probable error of a tabular 
value. In other words, such interpolated values are more precise than tab. 
values. The following table gives the values of r corresponding to several 
values of t. 

For t = or 1 r = ± 0.25 
" t = -At " A r = ± 0.23 
" t = •& « •& r = ± 0.20 
" * = A " A r = ± 0.17 

" * = A " A ' = ± °- 16 
" t = & " r = ± 0.15 
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As a statistical test of the accuracy of the laws of error expressed by eq's 
(8), (9), (10), 500 actual errors, e' = (1 — t)e 1 -\-ts 2 , corresponding to inter- 
polated values derived from a 5- place table of log's have been computed by 
means of a 7-place table. In each case t was taken equal to \, or the errors 
of interpolated mid-values were computed. The arguments used were the 
following numbers: 20005, 20035, 20065, 20105, 20135, 20165, 20205, 
etc., in the same order to 36635. Now for t — |, a = b = \, and hence by 
(8) and (10) 

f(e) = 4e-|-2 for values of e bet'n — J and 0, 
<p(e) = — 4e+2 " " " e " + | " 0. 

Prom the integral /(2±4s)de, the percentages of the whole number of er- 
rors to be expected between the limits and 0.1, 0.1 and 0.2, etc., up to 0.4 
and 0.5 may be computed. The following table gives a comparison of these 
percentages with those derived from the 500 actual errors. 



Between Limits. 


Percent of whole number of errors. 


By Theooy. 


By Trial. 


0.0 and 0.1 
0.1 " 0.2 
0.2 " 0.3 
0.3 " 0.4 
0.4 " 0.5 

0.0 " 0.15 


36 
28 
20 
12 

4 

50 


33.2 
30.2 
19.0 
13.2 
4.4 

51.4 



The average of the squares of all probable errors given by eq's (13) is 

(J)* | 2 Pfdt +4 f\l—i/WT^t)Jdt | = 0.594(|) 2 . 

2. The consideration of the probable errors of interpolated values of the 
second species presents some peculiar features, which are, perhaps, more in- 
teresting than important. These features arise from the fact already stated 
that e 3 in the formula £" = (1 — t)e 1 -}-te 2 -}-e z , is discontinuous for any val- 
ue of t. Thus, when t = |, e 3 must be either or J according as the tab- 
ulur difference is odd or even, and these two values of e 3 are equally likely 
to occur. For t = J therefore we may imagine any large number of result- 
ant errors, s", to be made up of two distinct systems of error, namely, 1st 
that which embraces all cases in which e 3 = 0, and 2nd that which embraces 
all cases in which e 3 = J. In the first system the probable error is ±0.15 
as given by the second of (13) for t = J, since e 3 = 0. In the second sys- 
tem all values of e" will be positive or negative according as e 3 is taken = 
-j-J or = —\; and the probable error will be + \ or — J according as e 3 = 
+Jor— \. 
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On the other hand, we may consider any large | 
number of values of t" , embracing both classes in I 
which £ 3 =0 and e 3 = \, as one system. In such 
a system, supposing the value for e 3 positive only, 
the actual errors will vary in magnitude between [ 
the limits — \ and -f 1. The law of facility for 
the errors in the system in which e 3 = only, is 
represented by the isosceles triangle ABO, the base 
AB being the axis of e and the altitude 00, being j 
the axis of f(e). Likewise, the isosceles triangle 
ODE represents the law of facility for the system I 
of errors in which s g = +J only, the origin being! 
at O with axes OE and 00 respectively. If the 
element-areas representing the probabilities of eq'l 
errors in the two systems be summed, the result j 
will be the probability area corresponding to the! 
combined system of errors. This area is repre- 
sented by the trapezoid AGDE. Therefore 

for values of s between — \ and 0, f{e) 
" " " e " « +hfh) 

" " " e " -H " +1, <p{e) 

If we consider all possible errors in this system with reference to magnil'de 
only, the probability-curve becomes a straight line joining the points whose 
coordinates are (0, 4) and (1, 0), such line being represented by EDF in the 
above Figure. The law of facility is 

<p(e) = 4(1 -«). 

To determine the probable error, r, we have the condition 

Ci(l—e)de = J C\(l~e)de = 1, whence 
J J 

4r — 2r a = 1 and 

r = 1 — 1^/2= +0.29. 
In order to verify by an actual test this result, which indeed differs from 
the result reached by a distinguished mathematician, see § 5, the actual er- 
rors of 500 interpolated mid-values from a 5-place table of logarithms have 
been computed by means of a 7-place table. These mid-values correspond 
to the arguments mentioned under §4, 1, above; and their actual errors cor- 
respond to e" in the formula s" = (1— t)s 1 -Ms 2 +e 3 , for t = J. Now by 
the law of error just derived, one-half of these actual errors ought to be less 
than 0.29. Likewise, from the integral /4(1 — e)d& 3 the percentages of the 
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whole number of errors to be expected between the limits and 0.2, 0.2 and 
0.4 etc., to 0.8 and 1 can be computed. The following table affords a com- 
parison of the theoretical and observed percentages. 



Between Limits. 


Per cent by Theo. 


Per cent by Obs. 


and 0.2 


36 


35.8 


" 0.29 


50 


49.8 


0.2 « 0.4 


28 


27.8 


0.4 " 0.6 


20 


18.6 


0.6 " 0.8 


12 


12.2 


0.8 " 1.0 


4 


5.6 



Again, suppose t — f . Then e 3 = 0, +£ or — J, and these values of e 3 
are equally likely to occur. Any large number of actual errors, t", for t = 
J may therefore be considered as a combination of three distinct systems 
corresponding respectively to the three different values of s 3 . For the sys- 
tem in which e s = only, <p (e) is given by equations (8), (9), (10) by put- 
ting a = § X J = f, and b = JX J = \. For the system in which e 8 = 
-f J only, f (e) will result from (8), (9), (10) by giving to a and b the values 
just mentioned, and substituting for e, £ + J; or in other words by moving 
the origin a distance + J to the left of its position for the system of errors 
in which e s = only. Likewise, for the system in which e 3 = — J only, 
f(e) will result by substituting for e in f(s) for the first system, e— J. See 
Fig., p. 149, t = £. If the areas representing the probabilities of equal er- 
rors in the three systems be summed, the probability-area for the original or 
combined system will result. This area will be a trapezoid whose base is 
J^-, whose side parallel to the base is f , and whose altitude is 3, the other 
two sides making angles with the base whose tangents are + f and — | re- 
spectively. Therefore 



for errors lying between 

It (l u li 



5- and -|, <p(e) = i(15+18e), 



— 7 
+* " 



+*, f(«) = 4(15-18e). 
Hence the probable error will result from 

f~5i(15-|-18e)de = f + %15— 18e)<fe = f, 

•/ — g- J +r 

which gives 

r = ±i[5 — 1/(12)] = ±0.26. 

Similarly it may be shown that for t = J, in which case e s = 0, -f \, 
J or — £, the possible errors being considered with respect to mag. only, 
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(p(e) = £(24 — 16e) for values of e between and \, 
f(e) = }(28 — 32e) " " " e " i " f, 
(pis) = J(16— 16e) " " " £ " | " 1, 
and r = ^[7-^(23)] = 0.28. 

Likewise for t = £, in which case £ 3 — 0, +i> +f> — $ or — i> i* mav 
be shown that 

(p(e) = 2£-(^-\-e) for values of £ between — ^j- and — ^, 

<p[S) = b £ — JQ- + T %, 

..(-A — 2 5/' 9 ff \ " " " ff " 4-1 « _L 9 

fl e J ¥ llT — V e +TT ±TTF> 

and r = ± r V[9— 1/(40)] = ±0.27. 

By extending the above process the probable errors in other and more 
complex cases might be determined. Without going into these details, 
however, it is evident that the probable error of an interpolated value of the 
second species will in general be greater than the probable error of a tabular 
value. For values of t such that 1-^i is large relatively to the largest tab- 
ular difference in the table used, £ s may become insignificant in comparison 
with (1 — ^-f-fej, and the probable error will be given in this case by for- 
mulas (13). At the limit, where I = 0, s s vanishes, and r = ±0.25. 

It is concluded, therefore, that the use of interpolated values of the sec'd 
species will not insure the highest precision possible with a given table. In 
precise computations, indeed, there can be no adequate reason for introdu- 
cing the error e 3 since the computer cannot avoid knowing its value. 

3. By reference to equation (5) it appears that the probable error of an 
interpolated value of the third species is the same as that of a tablar value. 
Hence there must be on the whole a loss of precision in using interpolated 
values of this species. 

4. For interpolated values of the fourth species the probable error will 
be given by equations (11) and (12) by making a = J and assigning to b 
the limiting value of s 4 , which is subject to the condition mentioned in §3. 
Whatever this limiting value of £ 4 may be, the probable error of the inter- 
polated value cannot be less than ± 0.25. When t is such that £ 4 may 
have any value between -f-f and — \ the probable error will be 

r = ±(1— -H/2) = ±0.29. 
At the limit, where t = and £ 4 = 0, 

r = ± 0.25. 

Historical Note. 

§ 5. Since at least three distinguished mathematicians have each consid- 
ered one or more of the problems discussed in this paper, without having 
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reached conclusions altogether accordant, a brief reference to their investi- 
gations may not be uninteresting. 

1. Bessel, 1838, in Nos. 358 and 359 of the Astronomische Nachrich- 
ten, in an article entitled Untersnchungen ueber die Wahrsoheinlichheit der 
Beobachtungsj 'elder, has investigated the laws of error discussed in §4, 1. 
This article is an admirable one and merits more attention than it seems to 
have received by writers on the theory of errors. In No. 529 of the Nach- 
richten, in a review of Captain Shortrede's 7-place tables of logarithms? 
Bessel has applied the laws of error defined by equations (8), (9), (10), to 
determine the probable error of an interpolated value of the first species. 
He observes that the proportional part may be computed by means of the 
tabular difference or by means of the true difference, giving rise to what we 
have called the first and third species. He then says 

Man bemerkt leicht, dass die erste methode die vortheilhaftere ist, aber 

ich erinnere niich nicht, dass der wahrscheinliehe Fehler ihres Resultats 

methodische untersucht worden wsere. 

Bessel does not raise the question of the independence of the errors of 

consecutive tabular values. His analysis leads to equations (13). He does 

not consider the probable errors of interpolated values of the second species. 

2. Bremiker, 1852, in the introduction to the Latin edition of his 6- 
place tables of logarithms has an important chapter entitled De Erroribus, 
Quibus Computationes Logarithmicae Afficiuntur. This chapt. aims to give 
a more comprehensive view of this subject than any work with which the 
writer is acquainted, and its omission in subsequent editions of Bremiker's 
tables would seem to have relegated it to unmerited obscurity. Among oth- 
er problems Bremiker considers that of the probable errors of interpolated 
values of the second species. (See his § 12.) His treatment of this problem 
is, however, quite erroneous owing to the fact that he has considered e 3 in 
s" = (1 — <)S] + t$ 2 -f- e 3 continuous between the limits -f- \ and — \ for 
any value of t. He says, p. 69, his f 1} f 2 , f 3 and e being the same as our 

e u s 2> £ 3 an< ^ t> respectively, 

Sumatur, f x et/ 2 esse errores logarithmorum tabulae L r et L 2 , ita in- 
telligendos, ut L x +/ x et _L 2 +/a sm * exacti logarithmorum valores, 
error logarithmi L 1 -\-sJ erit 

/i + </ 2 -/i) +fz 
si f 3 designat errorem, qui e producto decurtato e A proficiscitur, / vero 
putatur quilibet valorum seque probabilium, qui inter — T( == i) e t +r 
jacent Itaque primum quaeramus, quos valores summa 

(l-e)A+^+/, 
accipere possit, si pro illis / omnes valores, que esse possunt, aeque 
probabiles et intra — | et -f ^ jacentes ponantur. 
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He then states that the limits of this sum are — 1 and +1, and proceeds 
to determine the probable errors of interpolated values corresponding to sev- 
eral values of s = t. His results are very different from those reached in 
§4, 2, inasmuch as they increase gradually from a minimum ± 0.26 for 
t = ^-, to a maximum ± 0.29 for t — 0. The second of these probable er- 
rors is curiously absurd, since the inference from it is that an interpolated 
value infinitely near a tabular value is less accurate than a tabular value.* 
It may be remarked that the probable error ±0.26 for t = \ does not cor- 
respond to the percentages of actual errors given in § 4, 2. 

3. F. G. "W. Struve, 1860, in Arc Du Meridian, Tome I, p. 94, has 
considered the problem of the probable error of an interpolated value of the 
first species. He gives no analysis, but says briefly and somewhat obscurely 
L'erreur probable d'un logarithme donne' imme'diatement dans les ta- 
bles est 0.25 de la septieme de'cimale; l'erreur probable d'un loga- 
rithme de'duit par 1' interpolation entre deux chiffres tabulaires n' est 
que |/f de l'autre erreur. 
This statement may be interpreted thus. If in equation (2), e x and e a 
be supposed to conform to the ordinary law of error viz., f (£)=(/i.-j-j/a)e -7 ' 2fa 
the probable error corresponding to e' will be 

j[(l_<)2 + ^ _ i[l_ 2t+1t 2 Y\ 
Now the square root of the average of the squares of probable errors giv- 
en by this formula is 

iL/1(l— zt+2t 2 )dty = Wi = ±0.204. 
Although e x and e 2 do not conform to the ordinary law of error, it is 
interesting to note that the above formula gives results which do not differ 
widely from the true results given by equations (13). Struve makes no ref- 
erence to the writings of Bessel or Bremiker on this subject. 

*[Because e s must necessarily represent plus or minus 0, 1, 2, 3 or 4 
units in the last place figure of the interpolated value 

v' — {l—t)v i + tv 2 , 
and because any one of these digits is equally likely to occur, whatever the 
value of v' may be, the actual value of £ 3 is therefore independent of the 
value of v', depending only on the terminal digit in the numerical expres- 
sion for v', and is therefore equally continuous with s-^ and s 2 . But the 
probable error of the interpolated value will necessarily be increased by the 
addition of e 3 , and as its maximum occurs at the greatest distance of the 
interpolated value from mid-value, the absurdity above alluded to is not 
apparent. — Ed.] 



